草庐IT

Pygame Draw绘图

全部标签

单细胞/空间组利用scanpy实现Seurat的splitby分页绘图

在使用Seurat时,经常需要对不同分类的样本在同一画布上进行可视化,可以非常方便地通过其DimPlot()函数的goupyby和spliby等参数实现,例如对照组和实验组,如下:以正常和肿瘤组织作为区别但是在python的环境中,scanpy的sc.pl.umap()并没有这么灵活的参数。所以需要通过循环解决问题,sc.pl.umap中的color参数类似于Seurat的groupby,但其groups参数完全没有Seurat的splitby强大。所以我们可以通过python的Matplotlib包的plt.subplots()函数,结合循环将分组内容一一绘制,再多个分组图合并在一起,就实现

单细胞/空间组利用scanpy实现Seurat的splitby分页绘图

在使用Seurat时,经常需要对不同分类的样本在同一画布上进行可视化,可以非常方便地通过其DimPlot()函数的goupyby和spliby等参数实现,例如对照组和实验组,如下:以正常和肿瘤组织作为区别但是在python的环境中,scanpy的sc.pl.umap()并没有这么灵活的参数。所以需要通过循环解决问题,sc.pl.umap中的color参数类似于Seurat的groupby,但其groups参数完全没有Seurat的splitby强大。所以我们可以通过python的Matplotlib包的plt.subplots()函数,结合循环将分组内容一一绘制,再多个分组图合并在一起,就实现

3. 关于ggplot绘图公式复盘的一些问题2

  ggplot提供的几种重要几何图形函数geom_*如下:geom_smooth对数据进行拟合并展示拟合的结果和标准误;这个函数的重要参数是method。当method=NULL时,基于observations的数量(数据点的个数)选择默认方法:当observations的数量小于1000时,默认geom_smooth(method="loess");否则默认geom_smooth(method="gam",formula=y~s(x,bs="cs"))(ps:也可以写作geom_smooth(method="gam",formula=y~s(x))).geom_smooth(method=

3. 关于ggplot绘图公式复盘的一些问题2

  ggplot提供的几种重要几何图形函数geom_*如下:geom_smooth对数据进行拟合并展示拟合的结果和标准误;这个函数的重要参数是method。当method=NULL时,基于observations的数量(数据点的个数)选择默认方法:当observations的数量小于1000时,默认geom_smooth(method="loess");否则默认geom_smooth(method="gam",formula=y~s(x,bs="cs"))(ps:也可以写作geom_smooth(method="gam",formula=y~s(x))).geom_smooth(method=

统计绘图 | 这个能不能作为【火山图】?

1.火山图?之前瞥到了小明的跟着Cell学作图:R语言ggplot2作图展示差异表达的基因,但没有细看,根据印象做了下面这个图。差异表达结果合并awk'NR==1{print}FNR>1{print}'*vs0.txt>allvs0.txt在R中读取library(tidyverse)df_deg作图,其中不显示padj大于0.5的基因;df_deg的contrast列:1vs0|2vs0|3vs0|4vs0|5vs0,不同的组的差异表达结果;df_deg的change列:Up|Down|Stable,可根据log2FoldChange和padj定义;my_gene是我想突出的基因构成的向量。

统计绘图 | 这个能不能作为【火山图】?

1.火山图?之前瞥到了小明的跟着Cell学作图:R语言ggplot2作图展示差异表达的基因,但没有细看,根据印象做了下面这个图。差异表达结果合并awk'NR==1{print}FNR>1{print}'*vs0.txt>allvs0.txt在R中读取library(tidyverse)df_deg作图,其中不显示padj大于0.5的基因;df_deg的contrast列:1vs0|2vs0|3vs0|4vs0|5vs0,不同的组的差异表达结果;df_deg的change列:Up|Down|Stable,可根据log2FoldChange和padj定义;my_gene是我想突出的基因构成的向量。

R绘图:群体表型分布和相关关系

考察一个群体的多个表型或者一个表型的多个重复,我们想展示其分布和他们之间的相关关系可以使用柱状图和散点图(如下图所示)。test.png这幅图主要有两部分组成,一个是对角线上的柱状图,使用柱状图展示了每一个表型重复的分布;另一个就是对角线下面的散点图,用散点图展示两两之间的相关关系,并且用不同颜色表示点的密度,在上面标注其相关性。下面我们将使用R语言完成这幅图。对于这幅图我们可以先分别绘制其中每一个部分,然后使用图片组合、拼接函数进行整合:分图绘制首先导入数据,数据格式如下,每一行代表一个样本,每一列代表一个重复:>datadatapheno16rep1pheno16rep2pheno17ph

R绘图:群体表型分布和相关关系

考察一个群体的多个表型或者一个表型的多个重复,我们想展示其分布和他们之间的相关关系可以使用柱状图和散点图(如下图所示)。test.png这幅图主要有两部分组成,一个是对角线上的柱状图,使用柱状图展示了每一个表型重复的分布;另一个就是对角线下面的散点图,用散点图展示两两之间的相关关系,并且用不同颜色表示点的密度,在上面标注其相关性。下面我们将使用R语言完成这幅图。对于这幅图我们可以先分别绘制其中每一个部分,然后使用图片组合、拼接函数进行整合:分图绘制首先导入数据,数据格式如下,每一行代表一个样本,每一列代表一个重复:>datadatapheno16rep1pheno16rep2pheno17ph

在线绘图|凌波微课云平台搞定Lefse分析

1. 背景介绍LEfse分析即LDAEffectSize分析,可以实现多个分组之间的比较,还进行分组比较的内部进行亚组比较分析,从而找到组间在丰度上有显著差异的物种(即biomaker),主要是通过非参数因子Kruskal-Wallis秩和检验来实现的。LEfSe2.0在线工具在1.0工具基础上,增加了LDA判别分析结果文件(仅含差异显著)。增加差异特征图,以及所有结果图的重绘功能。工具链接:http://www.cloud.biomicroclass.com/CloudPlatform/SoftPage/LF22. 操作方法操作界面lefse分析需要输入的文件包含OTU丰度表和分组文件,OT

在线绘图|凌波微课云平台搞定Lefse分析

1. 背景介绍LEfse分析即LDAEffectSize分析,可以实现多个分组之间的比较,还进行分组比较的内部进行亚组比较分析,从而找到组间在丰度上有显著差异的物种(即biomaker),主要是通过非参数因子Kruskal-Wallis秩和检验来实现的。LEfSe2.0在线工具在1.0工具基础上,增加了LDA判别分析结果文件(仅含差异显著)。增加差异特征图,以及所有结果图的重绘功能。工具链接:http://www.cloud.biomicroclass.com/CloudPlatform/SoftPage/LF22. 操作方法操作界面lefse分析需要输入的文件包含OTU丰度表和分组文件,OT